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Abstract 

Nonperturbative Wilson coefficients associated with the leading chiral- 
symmetry-breaking operators in the operator product expansion of the pseu¬ 
doscalar QCD correlation function are derived. Implementation of the new, 
instanton-induced operators enables the corresponding spectral sum rule to 
reproduce the small pion mass scale, thereby reconciling it with Goldstone’s 
theorem. The same operators suppress the contributions of pionic resonances. 
Several predictions and structural insights from the new sum rule are dis¬ 
cussed. 
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The arguably most characteristic pion property is its mass well below all other hadronic 
mass scales. This special feature has long been understood as a consequence of spontaneous 
chiral symmetry breaking (SCSB) in the QCD vacuum [IJ, which renders the pion a (quasi-) 
Goldstone boson and manifests itself in the finite vacuum expectation values of chirality- 
changing operators like the quark condensate ( qq ). One would therefore expect the QCD 
sum-rule approach (2J, which explicitly links hadron properties by means of the operator 
product expansion (OPE) to such condensates, to provide a privileged source of information 
on the Goldstone nature of the pion and, in particular, to predict its small mass scale as a 
consequence of chiral-symmetry-breaking (CSB) condensates. 

Surprisingly, this expectation is not borne out by the existing analyses of the pseudoscalar 
sum rule. They fail to predict a pion mass below standard hadronic scales [|[, and the 
contributions from chirality-changing condensates are too strongly suppressed (by factors 
of the light quark masses) to have a notable impact on the resulting pion properties. This 
confronts us with a fundamental puzzle which threatens to unsettle the conceptual basis of 
the QCD-sum-rule approach: why is the conventional OPE of the pseudoscalar correlator 
practically blind to spontaneous chiral symmetry breaking, i.e. to exactly that element of 
QCD vacuum physics which profoundly affects the Goldstone boson channel? 

Additional contributions - neglected in the standard OPE but potentially enhanced in 
the pseudoscalar channel - could provide an attractive resolution of this puzzle. One such 
contribution, a low-dimensional power correction conjectured to originate from ultraviolet- 
sensitive physics, has been proposed in 0. Besides probably being small 0, however, 
this term is chirally invariant and therefore unlikely to significantly improve the pion mass 
prediction. Hard (so-called direct) instanton corrections to the unit-operator coefficient 
B-0. although of subtantial size, also proved unable to generate the physical pion mass scale 
||. The reason is probably the same as above: these contributions are chirally invariant (as 
their perturbative counterparts) and multiply the likewise chirally invariant unit operator. 

The above examples illustrate that information related to SCSB can enter the OPE 
only through chiral-symmetry-breaking condensates (since it originates from soft vacuum 
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physics). In search for missing physics related to such condensates, we will identify and 
calculate in this letter the leading nonperturbative corrections to the Wilson coefficients 
associated with the lowest-dimensional CSB operators. These new contributions arise from 
small instantons and have the potential to resolve the above-mentioned puzzle since in- 
stantons couple particularly strongly to pseudoscalar interpolating fields ]5|U and generate 
Wilson coefficients which are not suppressed by light-quark mass factors (in contrast to their 
perturbative counterparts). 

We start from the 7 r° correlation function 


fl (x) = (0\T j n o (x)j 7V 0 (0) |0), 


( 1 ) 


based on the pseudoscalar interpolator 

3ifi = J ^ {ui~f 5 u - di~f 5 d) , (2) 

which has an instanton-improved OPE (IOPE) (for a review see 0) of the general form 
n (/OPE) W 2) = j J d i xe m <0 |Tj„o (x)j%a (0) | 0 )<' opis > 

= E c 'o,(e 2 ;d(c , »w) (3) 


(Q 2 = — q 2 > 1 GeV, 11 < 1 GeV is the operator renomalization scale). 

The perturbative parts C of the Wilson coefficients, to O (a s ) for the unit operator 

G'n 

and to leading order for all remaining operators up to mass dimension d = 6 , are @,|l 0 |] 
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where 0 4j5 are the four-quark operators 

0 4 = —a {ua^ v t a u) [da^ u t a d) , 
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( t a = A“/2, where A a are the Gell-Mann matrices). 

Each term in the IOPE (^[) factorizes into contributions from hard modes with momenta 
|A:| > p, contained in the Cq , and from soft modes with \k\ < p in the operators O n . 
Despite widespread belief (based on asymptotic freedom) there are hadron channels where 
this does not even approximately amount to a factorization of perturbative and nonpertur- 
bative physics (at p ~ 1 GeV). Indeed, early conjectures f| and substantial recent evidence 
0111111 corroborate that (semi-) hard nonperturbative contributions due to small instantons 
are quantitatively relevant or even dominant in several hadronic correlators. 

The pseudoscalar correlator, in particular, is known to receive exceptionally strong direct- 
instanton contributions to the unit-operator coefficient, 


C[ W) (Q 2 ) = I dpn (ft) ^ 


™l,2 (p) + m h ( P ) 


(QpfKUQp) 


(9) 


ml, 2 (p) (p) 

(/, (/) refers to the (anti-) instanton closest to x, n (p) is the vacuum distribution of instan¬ 
tons with size p, and K\ (z) is a McDonald function ||T2[), which arise from the propagation 
of both quark and antiquark (ejected by (2j)) in the zero mode of the instanton held)] and 
can be obtained by means of semiclassical techniques [|.@. Due to interactions with am¬ 
bient, long-wavelength vacuum fields (including other instantons) [[L4| the quarks aquire an 
effective mass fh q% 2 (p) (the index indicates that two quarks are propagating in zero modes). 
The quantitative sum-rule analysis below only requires the value of m 9i 2 at the average in¬ 
stanton size p ~ 0.33 fm for which we adopt the recent estimate m 9i 2 (p) = fn q , 2 — 85 MeV 


obtained from instanton-liquid model (ILM) simulations of the pseudoscalar correlator |]T5 
The m 9i2 -dependence of the results will be discussed in 


All so far considered contributions to the IOPE coefficients are either associated with 
chirally-invariant operators or too strongly suppressed (note the factor m q in C ^ q ert ^) to 
generate more than minute corrections to the pion mass. Hence at this stage - which 


+ ^ (Q 2 ) h as the same momentum dependence as several instanton contributions (which arise 


from the analogous diquark loop) to baryon sum rules [11,13] 
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represents the current state of the art - the IOPE does contain virtually no information 
on the soft vacuum fields which are responsible for the spontaneous breakdown of chiral 
symmetry. Therefore it is not surprising that the corresponding sum rule is unable to 
generate the low mass scale which characterizes the Goldstone pion [[|. 

As we have argued above, there are reasons to believe that the missing information on 
SCSB is activated by nonperturbative contributions to the Wilson coefficients of chirally 
noninvariant operators which so far went unnoticed. Furthermore, direct instantons are 
promising candidates for such contributions since (i) they provide the leading nonpertur¬ 
bative deviations from asymptotic freedom, (ii) they are likely to play an important role 
in the dynamics of SCSB and in the strong flavor-mixing among pseudoscalar mesons ||, 
(iii) their small average size p < p~ 1 allows them to contribute strongly to the Wilson co¬ 
efficients, (iv) light-quark-mass suppression factors are absent, (v) the sensitivity of spin-0 
meson channels to instanton-induced short-distance physics is enhanced^] [|17]], and finally 
(vi) recent lattice measurements find the pseudoscalar correlator dominated by contributions 
from instanton-induced quark (quasi-) zero modes 

We are thus led to derive the instanton contributions to the Wilson coefficients associ¬ 
ated with the dominant (i.e. lowest-dimensional) chiral-symmetry-breaking operators of the 
IOPE, qq and gqaGq. They can be calculated as the leading terms in the semiclassical ex¬ 
pansion of the correlator around the (anti-) instanton in the background of long-wavelength 
quark and gluon vacuum fields. We postpone a detailed description of this calculation to 
0 and present here just the results, 


4rvc4/<w,)^ 
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(where iiq (a; b] z) is the confluent hypergeometric function JT2[) for the quark condensate 
coefficient and 


2 The same holds for spin-0 glueballs [18 
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for the coefficient associated with the mixed quark-gluon operator g s qaGq. This operator 
appears for the first time in the OPE of the pseudoscalar correlator. As expected, the 
above coefficients are not suppressed by small quark masses and can be several orders of 
magnitude larger than their perturbative counterparts. Since they arise from only one 
quark propagating in the zero-mode state (while the other, soft one contributes to the 
accompanying operator) we have denoted the corresponding effective mass fii q \ (p). 

Note that m 9j i does not equal fh q $ although both emerge from a mean-held picture of 
quark interactions with soft vacuum fields. The first main difference between the two is 
rooted in the fact that they arise from averages over the ensemble of vacuum fields (approx¬ 
imated, e.g., by instantons in the ILM). Due to the fluctuations in this ensemble one should 
not expect averages over more than one zero-mode propagator to factorize into separate 
averages over each propagator, and as a consequence m q ^ p fh q ^- This can be verified 
explicitly in the ILM framework [15| . 


The second difference is specific to the IOPE: while in the two-zero-mode contribution 
(j9|) the external momentum Q is shared between both quark lines, the full (i.e. maximal) 
Q hows through the one zero-mode quark line in (jlO]) and ©■ Now, in more complete 
treatments of the interactions with the vacuum background helds the effective masses will 
become momentum-dependent quark self-energies fh q (k). For momenta much larger than 
the chiral symmetry breaking scale, k A csb, these self-energies become insensitive to the 
soft CSB vacuum modes and approach the current quark mass, fh q (k —► oo) —> m q . Since 
the effective mean-held masses can be considered as momentum averages of such quark 
self-energies (or rfi q ^ = rh q (Qi) with Qi the typical momentum scale), one expects a scale 
hierarchy 


m q (Q = 0) > m q<2 > m qd > rn q 


( 12 ) 
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for all p. In order to get an idea of the size of m, h \ = m, q j (p, Q i), we will adopt the value of 
the quark self-energy M ( k ) in the large- N c approximation to the ILM [^0| at the momentum 
transfer Qi = 1.5 GeV (the mean value of the interval Q G [1, 2] GeV relevant for the sum 
rule below )[], 


% ~ M (Q = 1.5 GeV) = (Q) ~ 20 MeV, 


47r 2 p : 


(13) 


where M (0) ~ 0.3 GeV and <p (k) is a combination of modified Bessel functions given in 
sg. Since m 9i i ~ 0.02 GeV is rather close to the lower end of the admissible region (|T2|), 


the size of the instanton-induced coefficients (|T0| ) and (|lTD will reach about one half of their 
upper bound. The impact of different choices for l on predictions and stability of the 


sum rule will be discussed in Hl6 


Having calculated the IOPE of the pseudoscalar correlator up to d = 6, we now turn to 
the associated QCD sum rule. It will be convenient to rewrite the different parts n (A) ( Q 2 ) 
(.X G {pert, / + /,...}) of the correlator by means of the dispersion relation as 

ImnW (—s) 


n (X) (r) = H r n (x) (Q 2 ) = ^B r J 


ds- 


Q 2 


= —J ds!mn (A) (— s) e" 


(14) 


In Eq. ([14]) we have already applied the obligatoy Borel transform B r which improves IOPE 
convergence, removes subtraction terms and emphasizes the ground-state contribution to 
the correlator [Q. 

A spectral sum rule is then obtained by equating the IOPE description n( /OPE ) (r) 
of the correlator (in the r region where it is reliable, see below) to a standard hadronic 
representation n (p/l,en ' ) (r) whose spectral function consists of pole and duality-continuum 
parts: 


Im H {phen) (-s; s 0 ) = Im n (po * e) (-s) + Im n (cont) (-s; s 0 ). (15) 


3 The analogous estimate for m q ,2 yields the value — 85 MeV of Ref. [jH| at Q 2 /h-csB ~ 0.85, 
which might explain why fh q p is significantly smaller than typical constituent quark masses. 
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The effective threshold s o delimits the duality interval of the continuum into which we 
include, besides the standard OPE part, the instanton contributions: 


Imn (cont) (-s;s 0 ) 


Imn (OPE) (-s) + imn( /+f ) (-s) 


9(s- s 0 ) • 


(16) 


The instanton part will play an important role in the ensuing sum-rule analysis. In the pole 
(i.e. resonance) contribution we allow, besides the pion, also for its first excitation n', 

Im n (poie) (s) = nXlS (s - ml) + n\ 1,5 (s - m/,) , (17) 

where m n and m n > are the masses of the pion and the n', and A*. = y/2f^K (/| expt| = 93 MeV) 
with K = ml/ (m u + mj, ). Including the ir' resonance explicitly enables us to predict it’s 
strength A^, from the sum-rule analysis. Thus we can directly determine the quantitative 
impact of the tt 1 and decide whether it dominates the i r contribution (as suggested in [§]) or 
whether it can be absorbed into the duality continuum (as in other QCD sum-rules). 

Subtracting the continuum contributions of Eq. (|l^) from the IOPE, separately for each 
operator O n , we can write the sum rule as 


77 (r; s 0 ) = n (/OPS) (r) - (r; s 0 ) 
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77 


{pert) (r; s 0 ) + 77^ + ^ (t; s 0 ) 
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= \le-< T + 


(18) 


where the pole contributions are isolated on the right-hand-side and where we have defined 


i r°° 

77^ (r; s 0 ) = n<f } (r) - - / ds Im n<* } (-s) 6 (s - s 0 ) e— 
i r s ° / \ 

= -J ds (0 n ) IniCj 1 (-s) e~ ST . (19) 

It remains to calculate the imaginary parts of the Wilson coefficients in the timelike 
region from the explicit expressions for the CA ; given above. For the perturbative Wilson 
coefficients we find 
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The instanton-induced contributions can be obtained in closed form from the imaginary 


parts of the coefficients (|9|), (;1Q) and (ITT]). The unit-operator coefficient 


gives 
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The imaginary parts of the coefficients associated with chiral-symmetry-breaking operators 
have a more complex structure. They can be expressed in terms of the integrals 

V J+3 Ji (y/spv) 


"i 


hi,.j (s) = / dp- 
J o 


W«) = (24) 


>/W ’ J J o ' y/W ’ 

where J* (z) and V (2) are Bessel and Neumann functions fl2|]. For the instanton contribution 
to the quark-condensate coefficient (ill]) we find 


Imci 7+/) 


\ _ 2 6 7T“ / J n (p) p 4 , w ^ 2 3 7T 3 
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and for the mixed condensate coefficient (0) 

, 2V 


Im C^ I+I ^ 
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Note that both (^5|) and (^6|) receive strong contributions from s = 0 which significantly 
reduce the sp-dependence of the CSB terms lZ qq (r; so) and TZ gq Gaq (r; so) (cf. Fig. lb below). 

For the quantitative analysis of the sum rule ( |T8| ) we fix the IOPE parameters at standard 
values: A qcd = 0.2 GeV, (qq) = —0.0156 GeV 3 , ( a s G 2 ) = 0.04 GeV 4 , ( g s qaGq) = m^(qq), 
■tUq = 0.8 GeV 2 , (o^ + O^j = (567ra s /27) (qq) 2 . As in previous IOPE sum rules, we ap¬ 
proximate the instanton distribution as n (p) = nS (p — p) withQ n = 0.5 frn -4 and p = 0.33 
fin. The standard RG improvement, finally, amounts to scaling ^ O n ^ > £ 2 ~ 7n (o n ) with 

Tl ,a a G 2 0, 'Pqq 1, r )g s qaGq 1/6 aild 


f CO = ~2 ln ( rA2 ) 


2\ -4/9 


1 - 


290 
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256 ln [— ln (rA 2 )] 


729 ln (rA 2 ) 729 ln (rA 2 ) 


(27) 


4 Phenomenological, ILM and lattice evidence for these scales is discussed in [17] 
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to replacing a s by the (two-loop) running coupling, and to substituting p 2 —> 1/r. 

The resulting r- and s 0 -dependence of TZ (r; s 0 ) is plotted in Fig. la. Note that for 
r > 1.4 GeV -2 and so > 2.5 GeV 2 , TZ (r; so) becomes practically s 0 -independent. This 
is a consequence of the exponential exp (— st) in the integrands of (|T9| ) which renders the 
dispersive integrals insensitive to their upper integration limit So whenever exp (— Sqt) < 
ICG 2 . Therefore, So cannot be determined by the sum rule in this (r, So)-region. Another 
conspicuous feature of TZ (r; So) is the valley at intermediate So and small r, which can 
be traced to the instanton continuum contribution of the unit-operator coefficient. In this 
(r, So)-region the sum rule does not match since the positive slopes in r-direction do not 
fit the decaying exponentials of the pole contributions. In the so-region around 2 GeV 2 
(i.e. around the values which one would expect from duality arguments), however, TZ (r; so) 
shows an extended “mountain ridge” with the slow decay in r which matches an exponential 
containing a rather small mass. Thus we have qualitatively identified the area in which the 
sum rule will optimize, and we have found that the pion contribution must be sizeable. Both 
observations will be confirmed by the quantitative analysis below. 

Figure lb exhibits the r- and So-dependence of the new, instanton-induced CSB contri¬ 
butions to TZ (r; So), i.e. of 

7J< csb) (t; So ) = (t; s„) + U^a}, (t; so) . (28) 

The essential message of this plot is that in the physically meaningful (r, So)-region where 
So < 4 GeV 2 or r > 0.4 GeV~ 2 , the CSB contributions are monotonically increasing with 
r. In fact, those are the only relevant contributions to the IOPE whose slope is positive. 
Moreover, the slope at small r becomes maximal in the so-region around 2 GeV 2 where 
the sum rule optimally matches. Thus these CSB contributions overcome much of the 
negative slope originating from the chirally invariant operators (with both perturbative and 
instanton-induced coefficients) and thereby lower the prediction for the pion mass. This 
qualitative result demonstrates how the CSB condensates, “activated” by direct instantons, 
can reconcile the pseudoscalar sum rule with Goldstone’s theorem. As a byproduct, they 
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enhance the strength of the pion contributions relative to those from higher resonances. 

The main task of the quantitative sum-rule analysis is to find those values of the hadronic 
parameters to be determined for which both sides of ( |T8| ) optimally match in the fidu¬ 
cial T-region. The latter is obtained from the standard requirements that (i) the highest- 
dimensional operators should contribute less than 10% to 1Z (r; So), that (ii) multi-instanton 
effects should be negligible, and that (iii) the continuum contributions should be limited to 
maximally 50% of the total TZ (r; So). Since one sum rule cannot reliably determine all five 
hadronic parameters, we have chosen to fix those best known from experiment, m n o = 135 
MeV and m,i = 1300 MeV |22|, and to determine the values of the remaining onesQ, i.e. so, 
A,, and A^. Both sides of the optimized sum rule are shown in Fig. 2. Their good match con¬ 
firms the consistency and stability of the sum rule. The resulting values for the couplings are 
A 2 = 0.078 GeV 4 and A 2 , = 0.032 GeV 4 , while the continuum threshold becomes s 0 = 1.84 
GeV 2 . With f n = 93 MeV and m n = 135 MeV this implies rn u + ro^ = \/2f 7T rri^ /A^ ~ 9.0 
MeV (at n ~ 1 GeV), within the estimated range (6.75 — 16.2 MeV) of Ref. p2|. 


Fig. 2 also shows the individual contributions ^i^g 2 +o 4 +<>, Tz[ + ~\ lZ^ q + ^, and TZ^qGaq 
to the left-hand side of the sum rule. The instanton contributions can be seen to dominate, 
and the positive slope of the CSB parts indeed compensates most of the negative slope 
brought in by the chirally invariant ones. As expected, Tz{ + contributes with (strongly) 
negative slope. This explains why the sum rule of Ref. ||, which took exclusively this direct- 
instanton contribution (without it’s continuum part) into account, could not be stabilized 
without a strong resonance in the 1 GeV region. 

Implementing the CSB operators (and their continuum contributions) increases the rel¬ 
ative strength of the pion pole to (A 2 /A 2 ,)exp (rrv — nr n ) ~ 8 — 70 in the fiducial Borel 


5 Alternative analysis strategies will be considered in [ 161. We have checked, in particular, that 
the sum rule predicts the small pion mass scale, m n — 140 MeV, if A^ is fixed at it’s average 
phenomenological value A,r ~ 0.27 GeV 2 . 
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domain. As a consequence, the pion dominates while the n' contributions become negligi¬ 
ble. In conjunction with the relatively small value of the continuum threshold (s 0 = 1-84 
GeV 2 ) this suggests that the excited-state contributions can be absorbed into the duality 
continuum (as in practically all other sum-rule channels). The pion dominance seems natural 
in view of the exceptionally large mass difference between ground state and first excitation 
in the pseudoscalar channel. Moreover, it is consistent with lattice simulations of mesonic 
point-to-point correlators which find the pseudoscalar correlator well described by just the 
ground state pole and the duality continuum |23| . 

In summary, we have introduced instanton-generated Wilson coefficients associated with 
the leading chiral-symmetry-breaking operators in the IOPE of the pseudoscalar correlator. 
The new contributions are fully nonperturbative (both at soft and hard momenta) and 
supply previously missing information about spontaneous chiral symmetry breaking which 
reconciles the associated pseudoscalar sum rule with Goldstone’s theorem. As a consequence, 
this sum rule becomes the first in its channel which is able to reproduce the light mass scale 
of the pion. This resolves the puzzle stated in the introduction. Moreover, the chirally-odd 
operators suppress the contributions from higher-lying resonances, which can therefore be 
subsumed into the dispersive continuum. Both effects, as well as the stability of the sum rule, 
are enhanced by the instanton-induced continuum contributions. Additional implications of 
the new operators, e.g. for the calculation of the light quark mass values on the basis of 
pseudoscalar sum rules, will be discussed elsewhere 0. 


12 




REFERENCES 


[1] H. Pagels, Phys. Rep. 36, 137 (1978). 

[2] M.A. Shifman, A. I. Vainsthein, and V.I. Zakharov, Nucl. Phys. B147, 385, 448 (1979). 

[3] K.G. Chetyrkin, S. Narison, and V.I. Zakharov, Nucl. Phys. B550, 353 (1999). 

[4] M. Gockeler et al., Nucl. Phys. Proc. Suppl. 94, 571 (2001). 

[5] B.V. Geshkenbein and B.L. Ioffe, Nucl. Phys. B166, 340 (1980). 

[6] V.A. Novikov, M.A. Shifman, A. I. Vainsthein, and V.I. Zakharov, Nucl. Phys. B191, 
301 (1981). 

[7] E.V. Shuryak, Nucl. Phys. B214, 237 (1983); H. Forkel and M. Nielsen, Phys. Lett. B 
345, 55 (1995); A.S. Dorokhov, S.V. Esaibegyan, N.I. Kochelcv, and N.G. Stefanis, J. 
Phys. G 23, 643 (1997). 

[8] T.G. Steele, J.C. Beckenridge, M. Benmerrouche, V. Elias, A.H. Fariborz, Nucl. Phys. 
A624, 517 (1997). 

[9] H. Forkel, Instantons, OPE and Hadron Structure, and references therein, to be pub¬ 
lished. 


[10] C. Becchi, S. Narison, E. de Rafael, and F.J. Yndurain, Z. Phys. C 8, 335 (1981). 

[11] H. Forkel and M.K. Banerjee, Phys. Rev. Lett. 71, 484 (1993); H. Forkel and M. Nielsen, 
Phys. Rev. D 55, 1471 (1997); H. Forkel, F.S. Navarra, and M. Nielsen, Braz. J. Phys. 
31, 71 (2001). 

[12] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, Dover, New York 
(1965). 


[13] M. Aw, M.K. Banerjee and H. Forkel, Phys. Lett. B 454, 147 (1999); H. Forkel, |hep- 
| ph/0009275. 


13 






[14] M.A. Shifman, A. I. Vainshtein and V.I. Zakharov, Nucl. Phys. B163, 46 (1980). 


[15] P. Faccioli and E.V. Shuryak, Phys. Rev. D 64, 114020 (2001). 

[16] H. Forkel, to be published. 


[17] T. Schafer and E.V. Shuryak, Rev. Mod. Phys. 70, 323 (1998); for an elementary intro¬ 
duction to instantons see ff. Forkel, |hcp-ph/0009136 . 


[18] ff. Forkel, Phys. Rev. D 64, 034015 (2001), and |hcp-ph/0103204 . 


[19] T.L. Ivanenko and J.W. Negele, Nucl. Phys. Proc. Suppl. 63, 504 (1998); T. DeGrand 
and A. Hasenfratz, Phys. Rev. D 64 (2001) 034512; T. DeGrand, Phys. Rev. D 64 
(2001) 094508. 


[20] D.l. Diakonov and V.Yu. Petrov, Nucl. Phys. B272, 457 (1986). 

[21] C.A. Dominguez and E. de Rafael, Ann. Phys. 174, 372 (1987). 

[22] D.E. Groom et ah, Eur. Phys. J. C 15, 1 (2000). 

[23] M.-C. Chu, J.M. Grandy, S. Huang, and J.W. Negele, Phys. Rev. D 48, 3340 (1993). 


14 




I. FIGURE CAPTIONS 


1. Fig. la: The theoretical side of the sum rule, 7£(r;s 0 ), with m q; i = 20 MeV and 
fhq^ = 85 MeV. 


2. Fig. lb: The instanton-induced contributions of the chiral-symmetry-breaking opera¬ 
tors, JZS CSB "> (r; so) = 7^ + (r; so) + ( r > s o), to the theoretical side of the sum 

rule. 

3. Fig. 2: The right-hand-side (full line) of the optimized sum rule is compared to 

the theoretical side 1Z(t]s 0 = 1.84 GeV) (dotted). The contributions to 1Z from the 
perturbative (7 Z( pert \ short-dashed) and instanton-induced (Tl[ + ^ dot-dashed, lZ ^ q + 
dot-dot-dashed, short-dotted, their sum dashed) Wilson coefficients are also 

plotted separately. 
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